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Abstract 



This article is concerned with the mathematical analysis of the Kohn-Sham and 
■ extended Kohn-Sham models, in the local density approximation (LDA) and gener- 

^ I alized gradient approximation (GGA) frameworks. After recalling the mathematical 

^ ,'■ derivation of the Kohn-Sham and extended Kohn-Sham LDA and GGA models from 

the Schrodinger equation, we prove that the extended Kohn-Sham LDA model has a 
solution for neutral and positively charged systems. We then prove a similar result 
^ I for the spin-unpolarized Kohn-Sham GGA model for two-electron systems, by means 

of a concentration-compactness argument. 



^ ■ 1 Introduction 

a^ 

^ ■ Density Functional Theory (DFT) is a powerful, widely used method for computing approx- 

O ■ imations of ground state electronic energies and densities in chemistry, materials science, 

biology and nanosciences. 

According to DFT [IHllTS], the electronic ground state energy and density of a given 
molecular system can be obtained by solving a minimization problem of the form 



X 



inf {f{p)+ [ pV, p>0, ^eH\R^), [ 



p = N 



where N is the number of electrons in the system, V the electrostatic potential generated 
by the nuclei, and F some functional of the electronic density p, the functional F being 
universal, in the sense that it does not depend on the molecular system under consideration. 
Unfortunately, no tractable expression for F is known, which could be used in numerical 
simulations. 

The groundbreaking contribution which turned DFT into a useful tool to perform cal- 
culations, is due to Kohn and Sham who introduced the local density approximation 
(LDA) to DFT. The resulting Kohn-Sham LDA model is still commonly used, in particular 
in solid state physics. Improvements of this model have then been proposed by many au- 
thors, giving rise to Kohn-Sham GGA models [El EH El EH] , GGA being the abbreviation of 
Generalized Gradient Approximation. While there is basically a unique Kohn-Sham LDA 
model, there are several Kohn-Sham GGA models, corresponding to different approxima- 
tions of the so-called exchange-correlation functional. A given GGA model will be known 
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to perform well for some classes of molecular system, and poorly for some other classes. 
In some cases, the best result will be obtained with LDA. It is to be noticed that each 
Kohn-Sham model exists in two versions: the standard version, with integer occupation 
numbers, and the extended version with "fractional" occupation numbers. As explained 
below, the former one originates from Levy-Lieb's (pure state) contruction of the density 
functional, while the latter is derived from Lieb's (mixed state) construction. 

To our knowledge, there are very few results on Kohn-Sham LDA and GGA models in 
the mathematical literature. In fact, we are only aware of a proof of existence of a minimizer 
for the standard Kohn-Sham LDA model by Le Bris [13]. In this contribution, we prove 
the existence of a minimizer for the extended Kohn-Sham LDA model, as well as for the 
two-electron standard and extended Kohn-Sham GGA models, under some conditions on 
the GGA exchange-correlation functional. 

Our article is organized as follows. First, we provide a detailed presentation of the 
various Kohn-Sham models, which, despite their importance in physics and chemistry [24j . 
are not very well known in the mathematical community. The mathematical foundations 
of DFT are recalled in section [2l and the derivation of the (standard and extended) Kohn- 
Sham LDA and GGA models is discussed in section [3l We state our main results in 
section HI and postpone the proofs until section [5l 

We restrict our mathematical analysis to closed-shell, spin-unpolarized models. All 
our results related to the LDA setting can be easily extended to open-shell, spin-polarized 
models (i.e. to the local spin-density approximation LSDA). Likewise, we only deal with 
all electron descriptions, but valence electron models with usual pseudo-potential approx- 
imations (norm conserving |29|, ultrasoft [30], PAW [3]) can be dealt with in a similar 
way. 



2 Mathematical foundations of DFT 

As mentioned previously, DFT aims at calculating electronic ground state energies and 
densities. Recall that the ground state electronic energy of a molecular system composed 
of M nuclei of charges zi, zm (^a: S N \ {0} in atomic units) and electrons is the 
bottom of the spectrum of the electronic hamiltonian 



N N M 

=1 k=l ' * l<i<j<N 



where rj and Rfc are the positions in M'^ of the i*^ electron and the k^^ nucleus respectively. 
The hamiltonian Hj^i acts on electronic wavefunctions ^'(ri, cJi; • • • ; fat, cjAr), (jj G E : = 
{IT)) li)} denoting the spin variable of the i^^ electron, the nuclear coordinates {^k}Kk<M 
playing the role of parameters. It is convenient to denote by M| := x {||), ||)} and 
Xj := (rj,(Ti). As electrons are fermions, electonic wavefunctions are antisymmetric with 
respect to the renumbering of electrons, i.e. 

^'(Xp(i),--- ,Xp(Ar)) = e(p)^'(xi, • • • ,XAr) 

where e{p) is the signature of the permutation p. Note that (in the absence of magnetic 
fields) H]\j^ is real-valued if ^ is real-valued. Our purpose being the calculation of the 
bottom of the spectrum of i?Ar, there is therefore no restriction in considering real- valued 
wavefunctions only. In other words, if at can be considered here as an operator on the real 
Hilbert space 

N 

Hn = l\L\^,), 
1=1 
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endowed with the inner product 

(^-l^-^^ = / ^-Cxi,--- ,X7v)^''(xi,--- ,Xiv)dxi---(ixjv, 



where 



/ /(x)dx:= V / fir,a)dr, 

Jul Jr3 



i_ 

and the corresponding norm || • Wt-i^ = {•\-)'^^- It is well-known that Hj\[ is a self-adjoint 
operator on TIn with form domain 



TV 



^1) 



i=l 

Denoting by 

M 

Z = ^Zk 

k=l 

the total nuclear charge of the system, it results from Zhislin's theorem that for neutral or 
positively charged systems {Z > N), H]\f has an infinite number of negative eigenvalues 
below the bottom of its essential spectrum. In particular, the electronic ground state 
energy Xi{Hn) is an eigenvalue of Hj\f, and more precisely the lowest one. 
In any case, i.e. whatever Z and N, we always have 

XiiH^) = inf {(^iFTvl-f ), G Qn, imuj, = 1} • (2) 

Note that it also holds 

Xi{Hn) = inf {Tr (F^r), T G ^(Wtv), Ran(r) C Q^, < T < 1, Tr (r) = 1} . (3) 

In the above expression, SiTiiy) is the vector space of bounded self-adjoint operators on 
T-In, and the condition < T < 1 stands for < (*|r|^) < ||^'||^^ for all ^' e Hn- Note 
that if is a bounded-from-below self-adjoint operator on some Hilbert space Ti, with 
form domain Q, and if D is a positive trace-class self-adjoint operator on Tl, Tr (HD) can 
always be defined in ]R+U{-Foo} as Tr (HD) = Tr {{H - a)^ D{H - a)^2) + aTr (D) where 
a is any real number such that H > a. 

Prom a physical viewpoint, ([2]) and ([3]) mean that the ground state energy can be 
computed either by minimizing over pure states (characterized by wavefunctions or by 
minimizing over mixed states (characterized by density operators F). 

With any A''-electron wavefunction ^' G TIn such that H^'Ht^jv — 1 be associated 
the electronic density 

p*(r)=iVy' / |^(r,cr;x2, • • • ;xAr)p dx2 • • • dxAT. 

Likewise, one can associate with any A^-electron density operator F G S{Ti.iy) such that 
< F < 1 and Tr (F) = 1, the electronic density 

pr(r) = ^ / r(r, cj; X2, • • • , xat; r, a; X2, • • • ; xat) dx2 • • • dxAf 

(here and below, we use the same notation for an operator and its Green kernel). 
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Let us denote by 

M 

fc=i ' 



the electrostatic potential generated by the nuclei, and by 

N 

i 



Hl = -W^.,+ E ur^r (4) 



i=l l<i<j<N ' J' 

It is easy to see that 

(^'I/^atI^') = (^|//^|^)+ / p^V and Tr (HnT) = Tr (Hj^T) + [ prV. 

Besides, it can be checked that 

TZn = {p I 3^ G Q^, II^-IIt^^ = 1, pvi> = p} 

= {p\3Te SiHN), Ran(r) C Qat, < T < 1, Tr (r) = 1, pr = p} 



|p>o|^gfHm3), / 



p = N\. 
It therefore follows that 

In = inf + PG7^^| (5) 



inf|FL(p)+ f pV, pen^ 



7v| (6) 

where Levy-Lieb's and Lieb's density functionals [lH [15] are respectively defined by 

Fll(p) = inf {(^lif^l^), ^ G Qat, = 1, p^=p} (7) 

Fl(p) = inf {Tr (Hj^T), F E 5(?^Jv), Ran(r) C Q^, 

0<r<l, TV(r) = l, pr=p}. (8) 

Note that the functionals Fll and Fl are independent of the nuclear potential V, i.e. they 
do not depend on the molecular system. They are therefore universal functionals of the 
density. It is also shown in [15] that Fl is the Legendre transform of the function V In- 
More precisely, expliciting the dependency of In on V, it holds 

Fl(p) = sup (iNiV) - [ pv,ve + l°°(m3) 

from which it follows in particular that Fl is convex on the convex set TZn (and can be 
extended to a convex functional on L^(R'^) n L^(R^)). 

Formulae jS]) and jH]) show that, in principle, it is possible to compute the electronic 
ground state energy (and the corresponding groud state density if it exists) by solving 
a minimization problem on T^tv- At this stage no approximation has been made. But, 
as neither Fll nor Fl can be easily evaluated for the real system of interest (A'' inter- 
acting electrons), approximations are needed to make of the density functional theory a 
practical tool for computing electronic ground states. Approximations rely on exact, or 
very accurate, evaluations of the density functional for reference systems "close" to the real 
system: 

• in Thomas-Fermi and related models, the reference system is an homogeneous elec- 
tron gas; 

• in Kohn-Sham models (by far the most commonly used), it is a system of A'^ non- 
interacting electrons. 
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3 Kohn-Sham models 



For a system of N non-interacting electrons, universal density functionals are obtained as 
explained in the previous section; it suffices to replace the interacting hamiltonian Hj^ of 
the physical system (formula ^) with the hamiltonian of the reference system 



^ 1 

H% = -^-Ar,. (9) 



i=l 



The analogue of the Levy-Lieb density functional jT]) then is the Kohn-Sham type kinetic 
energy functional 

fKs{p)=mf{{^\H^^\^), ^ € Qn, W^W-Hr, = 1, Pvp =/9}, (10) 

while the analogue of the Lieb functional ((HI) is the Janak kinetic energy functional 

Tj{p) = inf {Tr {H%T), F G cS(H;v), Ran(r) C Q^, < T < 1, Tr (T) = 1, = p} . 

Let r be in the above minimization set. The energy Tr (H'^T) can be rewritten as a 
function of the one-electron reduced density operator Tr associated with T. Recall that 
Tr is the self-adjoint operator on L^(R|) with kernel 

Tr(x,x)=iV / r(x,X2,--- ,X7v;x',X2,--- ,XAr)(ix2---(ixAr. 

i(R|)JV-i 

Indeed, a simple calculation yields Tr (if^T) = Tr (— ^ArTr), where Ar is the Laplace 
operator on L^(M|) - acting on the space coordinate r. Besides, it is known (see e.g. [5]) 
that 

{T I 3r G SiTiN), Ran(r) C Qtv, < F < 1, TV (T) = 1, Tr = T, pr = p} 

= {T e 5(l2(m|)), < T < 1, Ran(T) C H\rI), Tr (T) = N, pr = p} , (11) 



where 



Hence, 



Pr{r) ■■= ^ T(r,cr;r,o-). 



Tjip) = inf |Tr ("^ArT^ , T G 5(L2(M|)), < T < 1, 

Ran(T) C H^M.^), Tr (T) = N, pr = pj. (12) 

It is to be noticed that no such simple expression for Tks(/o) is available because one 
lacks an A^-representation result similar to ([TT]) for pure state one-particle reduced density 
operators. In the standard Kohn-Sham model, Tks(/o) is replaced with the Kohn-Sham 
kinetic energy functional 

is a Slater determinant, p^ = p^ , (13) 
where we recall that a Slater determinant is a wavefunction ^ of the form 

^'(xi, . . . ,xjv) = ^=det((/)i(xj)) with (/>i G L2(M|), / </.i(x)0j(x) dx = Sij. 
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It is then easy to check that 



rKs(p) = inf(^V / |V</.i(x)|2dx, cD = (</,i,... ,</,^)G Wiv, = p] , (14) 



where we have set 



WN=f^^ = (.^i,--- An) I hHrI,), J^^ <^i(x),^,-(x) dx = 6ij^ 

and 

AT 

Note that for an arbitrary p E 7^ at, it holds 

Tjip) < TKsip) < Tks{p)- 

It is not difficult to check that p2|l always has a minimizer. If one of the minimizers T of 
(fT2]) is of rank N, then T = YliLi\4'i){4>i\ with $ = (0i,--- ,</'Af) G ^N, $ being then a 
minimizer of (fT3l) and rKs(p) = Tj{p). Otherwise, rKs(/o) > Tj{p). 

The density functionals Tks and Tj associated with the non interacting hamiltonian 
Hq are expected to provide acceptable approximations of the kinetic energy of the real 
(interacting) system. Likewise, the Coulomb energy 



representing the electrostatic energy of a classical charge distribution of density p is a 
reasonable guess for the electronic interaction energy in a system of electrons of density 
p. The errors on both the kinetic energy and the electrostatic interaction are put together 
in the exchange- correlation energy defined as the difference 

EM = Fll{p) - Tks{p) - J{p), (15) 



or 



EM = FUP)-Tj{p)-J{p), (16) 



depending on the choices for the interacting and non-interacting density functionals. We 
finally end up with the so-called Kohn-Sham and extended Kohn-Sham models 



tKS 



and 



7-EKS 

-'at 



= i^fj^y] / |V0,(x)|2(ix+ / p^V + Jip^) + E^,{p^), 

$ = ((/.I,--- ,(A7v) G Wjvj, (17) 

inf |Tr (-^ArT^ + j^^ prV + J{pr) + E^Pr), 

T G 5(L2(m|)), < T < 1, Tr (T) = iV, Tr (-ArT) < oo j, (18) 
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the condition on Tr (— ApT) ensuring that each term of the energy functional is well- 
defined. 

Up to now, no approximation has been made, in such a way that for the exact exchange- 
correlation functionals ( (fTSl) or (fTBI) ). I^^ = I^^^ = Xi{Hn) for all molecular system 
containing N electrons. Unfortunately, there is no tractable expression of E^dp) that can 
be used in numerical simulations. Before proceeding further, and for the sake of simplicity, 
we will restrict ourselves to closed-shell, spin-unpolarized, systems. This means that we 
will only consider molecular systems with an even number of electrons = 2Np, where 
Np is the number of electron pairs in the system, and that we will assume that electrons 
"go by pairs". In the Kohn-Sham formalism, this means that the set of admissible states 
reduces to 

^ = {(pia,(piP,- ■ ■ ,ifNpa,(fNpP) \ (Pi ^ H^O^^), / ipi(pj=Sij 

where a(|t)) = 1, a(|i)) = 0, /3(|T)) = and = 1, yielding the spin-unpolarized (or 

closed-shell, or restricted) Kohn-Sham model 



7-RKS 



Np 

^ = ((/.I,--- G (/^l(M^))^^ / ^iipj = 5ij, p^ = 2V|(^i|2l (19) 

where the factor 2 in the definition of p,p accounts for the spin. Likewise, the constraints on 
the one-electron reduced density operators originating from the closed-shell approximation 
read: 

T(r, IT), r', IT)) = T(r, |i), r', |i)) and T(r, |T), r', |i)) = T(r, |i), r', |T)) = 0. 

Introducing 7(r,r') = T(r, | T), r', | T)) and denoting by /0^(r) = 27(r, r), we obtain the 
spin-unpolarized extended Kohn-Sham model 

/^^^' = {^(7), l^lCNp} 

where 

f(7) = Tr (-A7) + / p^V + J{p^) + EUp^), 

and 

ICn, = {7 G cS(L2(m3)) I < 7 < 1, Tr (7) = Np, Tr (-A7) < 00} . 
Note that any 7 G /Cat^ is of the form 

-I- 00 

7 = ^ni\(t)i){(l3i\ 



i=l 



with 



i=l 1=1 



„ +00 +00 

(t)ieH^{R^), I (t)i(t)j=6ij, n, G[0,1], Vni = 7Vp, V ni||V</)i||i2 < 00. 
^R3 ~T ~T 

In particular, 

-I- 00 
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Let us also remark that problem (fT9]l can be recast in terms of density operators as follows 

I^^^ = {£{j), 7G/C^J (20) 

where 

= {7 e 5(L2(m3)) I 72 = 7, Tr (7) = Np, TV (-A7) < 00} 

is a the set of finite energy rank-A'p orthogonal projectors (note that fC^^ is the convex 
hull of Vnp)- The connection between (fTOll and (f20l) is given by the correspondence 

7 = X] 

1=1 

i.e. 7 is the orthogonal projector on the vector space spanned by the (pi. Indeed, as 
|V| = (-A)5, it holds 

Np Np Np 

TV (-A7) = Tr (|V|7|V|) = IW^Mh = E W^^^^Wh = E / l^-^^l'' 

i=l i=l i=l -^K^ 

Let us now address the issue of constructing relevant approximations for Exdp)- In 
their celebrated 1964 article, Kohn and Sham proposed to use an approximate exchange- 
correlation functional of the form 

Exc{p) = / 5'(p(r))dr (LDA exchange-correlation functional) (21) 

where p~^g{p) is the exchange-correlation density for a uniform electron gas with density 
p, yielding the so-called local density approximation (LDA). In practical calculations, it 
is made use of approximations of the function p ^ g{p) (from M+ to M) obtained by 
interpolating asymptotic formulae for the low and high density regimes (see e.g. [6]) and 
accurate quantum Monte Carlo evaluations of g{p) for a small number of values of p [1]. 
Several interpolation formulae are available [23l [22j [3l] , which provide similar results. In 
the 80's, refined approximations of -Exc have been constructed, which take into account the 
inhomogeneity of the electronic density in real molecular systems. Generalized gradient 
approximations (GGA) of the exchange-correlation functional are of the form 

Exc{p) = I h{p{r), -|V-y/ /j(r)|2) dx (GGA exchange-correlation functional). (22) 

Contrarily to the situation encountered for LDA, the function (p, k) i-^ g{p, k) (from M+ x 
M+ to M) does not have a univoque definition. Several GGA functionals have been proposed 
and new ones come up periodically. 



Remark 1. We have chosen the form / fUj) for the GGA exchange-correlation functional be- 
cause it is well suited for the study of spin-unpolarized two electron systems (see TheoremlE 
below). In the Physics literature, spin-unpolarized LDA and GGA exchange-correlation 
functionals are rather written as follows 



EM = EAp) + Ec{p) 



with 



EAp) = [ p(r)ex(p(r))Fx(sp(r))dr (23) 
Ecip) = I p{r)[e,{r,{r))+H{r,{r),tp{r))]dr. (24) 
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In the above decomposition, is the exchange energy, is the correlation energy, and 
Cc are respectively the exchange and correlation energy densities of the homogeneous electron 
qas, r^fr) = (Ivrpfr)) ^ is the Wiqner-Seitz radius, So(r) = ^ — ^IZ^Wi 

^ ^3 ^ ' ' 2(37r2)3 p(r)i ' 

dimensional) reduced density gradient, tp{r) = — ^ i) ^ ^'^P^^^^ ig ffig correlation gradient. 

Fx is the so-called exchange enhancement factor, and H is the gradient contribution to the 
correlation energy. While ex has a simple analytical expression, namely 

ec has to be approximated (as explained above for the function g). For LDA, F^ is every- 
where equal to one and H = 0. A popular GGA exchange- correlation energy is the PBE 
functional [2U], for which 

FJs) = 1+ 



1 + fiu ^s"^ 

the values of the parameters fj. ~ 0.21951, u ~ 0.804, 6 = 7r~^(l — In 2) and v = 37r~^/i 
following from theoretical arguments. 

4 Main results 

Let us first set up and comment on the conditions on the LDA and GGA exchange- 
correlation functionals under which our results hold true: 

• the function g in (f2T]) is a function from M+ to M, twice differentiable and such 
that 

5(0) = (25) 

g' <0 (26) 

30 < /3_ < /3+ < ^ s.t. sup j^y' < oo (27) 
31<a<^ s.t. limsup^<0; (28) 

• the function h in (|2T]) is a function from Mj_ x M+ to M, twice differentiable with 
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respect to the second variable, and such that 

h{0, n) = 0, Vk G M+ (29) 

dh , , 

a? £ (30) 



2 

30 </?_</?+< - s.t. sup 



dh 



< oo (31) 



^3 h(p,K) , , 

31 < a < - s.t. hmsup ^ < (32) 



2 



(p,k)~*(o+,o+) P' 



dh 

30<a<6<oo s.t. V(p,k) G M+ X R+, a < 1 + — (p, k) < 6 (33) 

OK 

dh d'^h 

V(p,k) GM+ xM+, 1 + —{p,k)+2k^{p,k)>Q. (34) 



Conditions (|25]) - (l28l) on the LDA exchange-correlation energy are not restrictive. They 

are obviously fulfilled by the LDA exchange functional (5x^^(p) = ~\ (f)^ P^)i ^'^d are 
also satisfied by all the approximate LDA correlation functionals currently used in practice 
(with a = I and /?_ = /3+ = |). We have checked numerically that assumptions ([29]) - ([34l) 
are satisfied by the PZ81 functional defined in |23| . 

Remark 2. Our results remain true if (2^) and are respectively replaced with the 
weaker conditions 

^/ ^ / 2 max(0,g^(p)) 

3- < /3_ < /3+ < - s.t. sup — J, — < oo 

-3 peK+ /r- + p/^+ 

and 

max fo, —(/),/«) 
3- </?:</?+<- s.l sup ^<oo. 

(p,K)gR+xR+ p^- + 

As usual in the mathematical study of molecular electronic structure models, we embed 
([20]) in the family of problems 



lA = inf{f(7), 7G/Ca} (35) 

parametrized by A G M+ where 

/Ca = {7 G S{L\R^)) I < 7 < 1, TV (7) = A, Tr (-A7) < 00} , 
and introduce the problem at infinity 

=inf{f°°(7), 7e/CA} (36) 

where 

fKS(7) = Tr (-A7) + J{pj) + E,,{p^). 
The following results hold true for both the LDA and GGA extended Kohn-Sham models. 



Lemma 1. Consider f35]) and (36\) with E^c given either by (21]) or by a2^) together with 
the conditions ^M-^M) or (M-IEE. Then 



10 



1. 1^ = 1^ = and for all X > 0, -oo < h < < 0; 

2. the functions X ^ I\ and X are continuous and decreasing; 

3. for all < fi< X, 

h<I, + ir^,. (37) 
Our main results are the following two theorems. 

Theorem 1 (Extended KS-LDA model). Assume that Z > N = 2Np (neutral or positively 
charged system) and that the function g satisfies \25\}- (28\} . Then the extended Kohn-Sham 
LDA model fg^) with i?xc given by has a minimizer 70. Besides, 70 satisfies the 
self- consistent field equation 

'yO = X{--oo,ey){Hp^^)+S (38) 

for some ep < 0, where 

Hpfo = + ^ + P^o * kr' + g'ip^o), 

where X{-oo,ep) is the characteristic function of the range (—00, ep) and where 6 G 5(-L^(R^)) 
is such that < 6 < 1 and Ran{5) = Ker{Hp^^ — ep). 

Theorem 2 (Extended KS-GGA model for two electron systems). Assume that Z > N = 
2Np = 2 (neutral or positively charged system with two electrons) and that the function 
h satisfies ^MM)-^3^. Then the extended Kohn-Sham GGA model (E^i with E^c given by 
^2^} has a minimizer 70. Besides, 70 = \4>){4>\ where (p is a minimizer of the standard 
spin-unpolarized Kohn-Sham problem f7g|) for Np = 1, hence satisfying the Euler equation 

_idiv (^(^l + |^(p^,|V</.|2)^ V^) + (y + /5<^*|r|-^ + |^(p^,|V0|2)^ = (39) 

for some e < 0, where = 20^. In addition, (j) G C'^'"(M^) for some < a < 1 and decays 
exponentially fast at infinity. Lastly, (j) can be chosen non-negative and (e, (j)) is the lowest 
eigenpair of the self-adjoint operator 



-l""' ((l + S(^'*.|V*l"))v-)+V + P**|r|-' + 3^ 

We have not been able to extend the results of Theorem [2] to the general case of Np 
electron pairs. This is mainly due to the fact that the Euler equations for ([35]) with E^c 
given by ([22]) do not have a simple structure for Np > 2. 

5 Proofs 

For clarity, we will use the following notation 
E}:?Hp)= [ 9{p{r))dT 

Eg''Hp)= I h{p{v)^-\V^p{v)\^)dr 
Jr3 2 

£^^'^{j)=Tt{-Aj)+ [ p^V + J{p^)+f g{p,{v))dv 

Jrs Jr3 

^GGA(^) ^ ^ /■ p^y^j^p^^^ f /j(p^(r),i|Vv7^(r)|2)dr. 

JR3 Jr3 ^ 

The notations E^dp) and f (7) will refer indifferently to the LDA or the GGA setting. 
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5.1 Preliminary results 



Most of the results of this section are elementary, but we provide them for the sake of 
completeness. Let us denote by Si the vector space of trace-class operators on L^(M^) 
(see e.g. [25]) and introduce the vector space 

W = {7 G ©1 I |V|7|V| G Si} 

endowed with the norm || • ||-^ = TV (| • |) + Tr (||V| • |V||), and the convex set 

/C = {7 G S{L^{R^)) I < 7 < 1, Tr (7) < 00, Tr (|V|7|V|) < 00} . 



Lemma 2. For all 7 G /C, -^/p^ G H (M ) and the following inequalities hold true 
1 



^llV^P^IIi^ <Tr(-A7) (40) 



< J(p^) < C(Tr 7)2(Tr (-A7))2 (41) 
-4Z(Tr 7)5(Tr {-A-f))^ < / p^V < (42) 

/ -, P- 3/3_ ,3+ 3/3+ \ 

-C f (TV 7)i--(Tr (-A7))— + (TV j)'~-{Tt (-A7)) — j < E^,{p^) < 0(43) 



^(7) > - (^(Tr (-A7))2 - 4Z(Tr 7)2 j - 8Z^Tr 7 

/ 2-,3_ 2-/3+ \ 

-C (Tr 7) 2-3."- + (TV 7) ^-■'^+ (44) 



2-l3_ 2-/3, 



£^{l) > ^Tr (-A7) - C (^(Tr 7) + (Tr 7) j , (45) 

for a positive constant C independent ofj. In particular, the minimizing sequences of |2 
and those of (3E\} are bounded in TL. 



Proof. Any 7 G /C can be diagonalized in an orthonormal basis of L^(M^) as follows 

+00 

7 = ^nj|(/))((/)i| 

i=l 

with n, G [0, 1], G H^{M?), J^, = 6ij, TV (7) = EtT < ^ and Tr (-A7) 
E^tT^^llV^.lli^ <oo. As 



ni(l)i\/(l)i 



i=l 



1=1 



iO]) is a straightforward consequence of Cauchy-Schwarz inequality. Using Hardy-Littlewood- 
Sobolev [16], interpolation, and Gagliardo-Nirenberg-Sobolev inequalities, we obtain 

J{Pj) < c,\\p,f e < c,\\p41\\p41, < C2\\p4l\\v^h^. 

Li ^ 
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Hence JH]), using (PO]) and the relation ||/0^||li = 2Tr (7). It follows from Cauchy-Schwarz 
and Hardy inequalities and from the above estimates that 




< 2||p^||1i||VVpVIIl2 < 4(TV 7)i(Tr (-A7))l. 



Hence g2]). Conditions (l25|)-(l28D for LDA and (I29|)-(l32|) for GGA imply that E^^ip) < 
and there exists 1 < p- < p-^- < ^ {p± = 1 + P±) and some constant C € M+ such that 

yp G fC, \E^,ip)\ <c( [ pP-+ [ /+V (46) 

from which we deduce ([43]l . using interpolation and Gagliardo-Nirenberg-Sobolev inequali- 
ties. Lastly, the estimates (|44|) and (l45l) are straightforward consequences of ([4T]) - ([43]) . □ 



Lemma 3. Let A > and 7 G ICx. There exists a sequence (7n)neN such that 
1. for all n G N, 7„ G K,\, 7„ is finite-rank and Ran{'^n) C C^(IR'^); 
■2. (7n)neN converges to 7 strongly in TL; 
^- (■^P7„)neN converges to ^Jp^ strongly in /f^(M^); 

4- (P7n)n6N \/P'yn)nen Converge almost everywhere to p^ and Vy/p^ respectively. 

In particular 

lim £{jn)=£{l) and lim £^ {jn) = S"^ il) ■ (47) 

n— >oo n—foo 

Proof. Let 7 G /Ca- It holds 

+00 

7 = 

i=l 

with ni G [0,1], G //i(IR^), /R3</'i</'i = Tr (7) = = X and Tr (-A7) = 

E^tT^.IIVc/'.lli, <oo. 

We first prove that 7 can be approached by a sequence of finite-rank operators. Let 
A'o G N such that < nAr^ < 1 (if no such A'o exists, then 7 is finite-rank and one can 
directly proceed to the second part of the proof). For all G N, we set 



IN 



N / ^ \ 

'^ni\(j)i){(pi\ + I A - I 

i=l \ i=l / 



I'^A'oX'AtVoI 



For N large enough, 'yjy G JCx, and the sequence (7Ar) obviously converges to 7 in TC. 
Besides, (/O^^) converges a.e. to p^ and 



N 



i=N+l 



i=l 



Hence the convergence of (p^^) to p-y in LP(]R ) for all 1 < p < 3. Besides, for all > A'o, 

2 



N / ^ \ 

i=l,i^No \ i=l / 



7jvI 



N 



N 



< 2J2n^\V(|)^\^+2X\V^No? 



i=l, ij^No 



+ njvo + A - ^ 



ni \(pNo\ 



i=l 



i=l 
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Using Lebesgue dominated convergence theorem, we obtain that the sequence (||Vy^p^||^2) 
converges to ||Vy^||^2, from which we deduce that {^/p^) converges to y/p^ strongly in 

The second part of the proof consists in approaching each cpi by a sequence of regular 
compactly supported functions. For each i, we consider a sequence ((/'i,fc)fegN of functions 
of C~(M3) such that 

• supp{<pi^k) C supp{4>i) and /^a (t>i,k(t>j,k = ^ij for all k, 

• (0j,fc)fcGN converges to 0j strongly in //^(M^) and almost everywhere, 

• there exists hi G L^(R'^) such that \V(j)i^k\ < hi for all k. 
It is then easy to check that the sequence {jN,k)k<^N defined by 

N / ^ \ 

lN,k = ^^ni\(l)i^k){4'i,k\ + I X^^i I \'pNo,k){(t)No,k\ 

i=l \ i=l / 

converges to in and is such that {^/P^^)keN converges to ^fp^ strongly in //^(M^). 

One can then extract from (7Ar,fc)(Ar,fc)gN*xN a subsequence (7n)neN which converges to 
7 in W and is such that (y^p^)neN converges to ^Jp^ strongly in i?^(M^), and there is no 
restriction in assuming that (p7„)neN and ^ ^fp^n&\ converge almost everywhere to 
and V ^/p^ respectively. 

The linear form 7 i— > Tr (— A7) being continuous on Ti and the functionals u ^ j^j. v?V 
and u ^ Jiu^) + Ey^du^) being continuous on i7^(]R^), (l47ll holds true. □ 

5.2 Proof of Lemma [U 

Obviously, Iq = 1^ = and h < for all A € M+. 

Let us first prove assertion 3. Let 0</i<A,e>0 and 7 G fC^ such that < <?(7) < 
+ e. It follows from Lemma [3] that there is no restriction in choosing 7 of the form 

TV 

7 = ^ni\(l)i){(j)i\ 

i=l 

with < nj < 1, X^i^i ^« — A*) {4'i\4'j) = ^ij and (j)i G C^(M^). Likewise, there exists 

N' 

1=1 

with < < 1, EZiK = a - /X, {(l)'i\(t>'j) = 5ij and G C^{R^), such that < 
< -^A^^ + e. Let e be a unit vector of M? and Tq the translation operator on L^(M^) 
defined by Tq/ = /(• - a) for all / G L'^{R^). For n G N, we define 

7n = 7 + Tnel'T-ne- 

It is easy to check that for n large enough, 7„ G /Ca and 

IX < £{ln) < £{l) + £^{l) + D{p^, T^,py) <I^ + + 3e, 

where 

D(p,/):= / / «^*<ir', 

Jr3 7]83 |r-r'| 
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Hence jS] 

Making use of similar arguments, it can also be proved that 

IT<I^ + ir-^.- (48) 

Let us now consider a function (j) € C^(M^) such that \\4'\\l'^ = 1- For all cr > and all 
< A < 1, the density operator 7(j,a with density matrix 

7^,A(r, r') = AfJ^ cpiar) 0(crr') 

is in /Ca- Using (f28ll for LDA and (f32l) for GGA, we obtain that there exists 1 < a < |, 
c > and ctq > such that for all < A < 1 and all < o" < cio, 



Therefore < for A positive and small enough. It follows from ((37l) and ([48]) that the 
functions X ^ Ix and A i-^ are decreasing, and that for all A > 0, 

-oo < /a < /~ < 0. 

To proceed further, we need the following lemma. 

Lemma 4. Let A > and (7n)nGN be a minimizing sequence for (E^. Then the sequence 
{Pjn)neN cannot vanish, which means that 

3R > s.t. lim sup / p^^ > 0. 
The same holds true for the minimizing sequences of /i36\) . 



Proof. Let (7n)neN be a minimizing sequence for (|35ll . By contradiction, assume that 



Vi? > 0, lim sup / pn = 0- 

"^°°a;gR3 Jx+Br 



Let 1< p < |. For p > such that ^ G H'^{M.^), it holds for all k G Z^, 

k+Bi \Jk+Bi J \Jk+Bi / \Jk+Bi / \J k+Bi 

(where the constant Cp does not depend on k). We therefore obtain 

/ / < E / / 

JR3 ^ Jk+Bi 

p-1 



keZ^ VJfc+Bi / \Jk+Bi / 

< ScJsnp f pY ( f {p+ f \V^\'] 



Hence, for all 7 G /C, 

p-i 



I p^<16Cpfsup / 



P7 ) II7IIH- 
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As we know that any minimizing sequence of (f35l) is bounded in H, we deduce from the 

5 

3' 



above inequahty that for all 1 < p < |, 



hm / pP^=0. 
In particular, it follows from ([46]l that 

lim / E^cip^J = 0- 

Let us now fixl<p<|, e>0 and R > such that \V\ < eX^^ on Bf^. For n large 
enough, we have 

/ P,.V < [ p,jV\+ [ p,jV\<([ iVlA^'d' + P,.<2e. 

JR^ J Br JbZ K-JBr J \JBr J ^ JB-j^ 



Therefore 



As, 



lim / p^„y = 0. 



we obtain that Ix > 0. This is in contradiction with the previously proved result stating 
that Ix < 0. Hence {P'y„)n€N cannot vanish. The case of problem ([361) is easier since the 
only non-positive term in the energy functional is £'xc(p)- □ 

We can now prove that Ix < Ix^- For this purpose let us consider a minimizing sequence 
(7n)ngN for ([36]) . We deduce from Lemma |4] that there exists rj > and R> 0, such that 
for n large enough, there exists x„ S such that 

/ P7n > V- 

Jxn + Bn 

Let us introduce 7n = Tx^-xr^lnTxr^-xi- Clearly 7„ G Kx and 



Thus, 



£{ln)<£°^{ln) ^ 



J ^ TOO ^ TOO 



It remains to prove that the functions X ^ Ix and A I^ are continuous. We will deal 
here with the former one, the same arguments applying to the latter one. The proof is 
based on the following lemma. 

Lemma 5. Let (afc)fcgN be a sequence of positive real numbers converging to 1, and {pk)ke'N 
a sequence of non-negative densities such that (y^)jtgp^ is bounded in H^{M.^). Then 

lim {E^ciakPk) - E^dPk)) = 0. 

k^oo 
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Proof. In the LDA setting, we deduce from ([28]l that there exists 1 < P- < p+ < | and 
C £ M+ such that for k large enough 

\E^F^{c.kPk) - E^^Hpk)\ < C |a, - 1| / {pI' +pI^). 

In the GGA setting, we obtain from (fSTI) and ((33]) that there exists 1 < p- < p+ < | and 
C S M4- such that for A; large enough 

\Eg''\a,p,) - E^,^Hpk)\<C\a,-l\ [ (p^ + + | V^P^p). 

As (y^)fegN is bounded in H^CK^), {pk)keN is bounded in LP(R^) for all 1 < p < 3 and 
{Vy/pk)keN is bounded in (L^(R^))^, hence the result. □ 

We can now complete the proof of Lemma [TJ 

Left- continuity of X ^ I\. Let A > 0, and {Xk)kGn be an increasing sequence of positive 
real numbers converging to A. Let e > and 7 G ICx such that 

h < f (7) <h + ^- 

For all A; € N, 7fe = AfcA~"^7 is in /C^fe so that 

VA; G N, Vn G N, h < h, < £{lk). 

Besides, 

^(7.) = ^Tr (-A7) + ^ J^^ p,V + ^ J(p,) + i?., (^p,) ^-^£:(7) 
in virtue of Lemma [3 Thus 

/a < /a, < /a + e 

for k large enough. 

Right-continuity of I\. Let A > 0, and (Afc)^^^ be an decreasing sequence of positive 
real numbers converging to A. For each A; G N, we choose 7/c G 1C\^ such that 

1 

k' 



h, < Silk) < h, + T- 



For all A: G N, we set 'jk = AA^, ""^7^. As ik G IC\, it holds 

h < £{lk) = ^Tr (-A7fc) + A /■ p^^v + ^ J(p^J + i^xc f ^P7. ) ■ 

As (7A:)fceN is bounded in TL and {^pj^)keN is bounded in i?^(M^), we deduce from LemmaE] 
that 

lim {£{%) - f (7fc)) = 0. 



Let e > and k^ > 2e such that for all k > k, 



mk)-nik)\<^. 

Then, 

VA: >k„ h-e< h, < h- 

This proves the right-continuity of A 1— > /a on R+\{0}. Lastly, it results from the estimates 
established in LemmsEJthat 

lim Ia = 0. 

A-*0+ 
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5.3 Proof of Theorem [U 

Let us first prove the following lemma. 

Lemma 6. Let (7ri,)neN o, sequence of elements of IC, bounded in 7i, which converges to 
7 for the weak-* topology ofTC. //lim„^ooTr (7„) = TV (7), then ip-y„)neN converges to p-y 
strongly in LP(M'^) for all I < p < 3 and 

£:LDA(^) < liminf£:LDA(7„) and S^^^'^^i-f) < liminf £:L^^'~(7„). 

n— >oo n— >oo 

Proof. The fact that (7n)jieN converges to 7 for the weak-* topology of Ti. means that for 
all compact operator K on L^(M^), 

lim TV (7„K) = TV (7K) and lim Tr (|V|7„| V|i<:) = TV (| VI7I Vl-fsT). 

For all W G C^{R'^), the operator (1 + |V|)~^VF(1 + |V|)-^ is compact (it is even in the 
Schatten class &p for all p > | in virtue of the Kato-Seiler-Simon inequality |27j), yielding 

/ py^W = 2Tr(7„Ty) = 2Tr((l + |V|)7„(l + |V|)(l + |V|)-iVF(l + |Vri) 

^ 2Tr ((1 + |V|)7(1 + |V|)(1 + \V\)-'^Wil + = 2Tr (7!^) = / p^W 



Hence, (p7„)ngN converges to pj in V'iW^). As by ([40]l . (■^P7„)neN is bounded in 
it follows that (y^p^)neN converges to y/p^ weakly in H^{M.^), and strongly in L|'^^(]R^) for 
all 2 < p < 6. In particular, {^yp■y„)n€N converges to y/p^ weakly in L^(M^). But we also 
know that 

lim \\^/p^\\l2 = lim / pj„ = 2 lim Tr (7„) = 2TV (7) = / p^ = H^H^a. 

Therefore, the convergence of {^/p^)nGN to y/p^ holds strongly in L^(M^). By an elemen- 
tary bootstrap argument exploiting the boundedness of (y^p^)neN in -ff^(M^), we obtain 
that {^yP'y„)nGN couvcrgcs strougly to y/p^ in LP(M^) for all 2 < p < 6, hence that (p7„)neN 
converges to p^ strongly in LP(M^) for all 1 < p < 3. This readily implies 

lim / P-y„V = pyV 
lim J{pjJ = Jip-r) 

n—*oo 

lim ^LDA( j^^LDA( 
n— >oo 

Lastly, for any orthonormal basis {•ipk)k€N* of L^(M^) such that xpk £ -ff^(M^) for all k, we 
have 

+00 

Tr(|V|7|V|) = j;(V'fc||V|7|V||V'fc) 
fc=i 
+00 



^Tr(7(||V|^fc)(|V|V'.|)) 

k=l 
+00 

^ lim Tr (7„(||V|V'fc)(|V|V'fc|)) 



fc=l 



< 



liminf j;TV(7„(||V|Vfc)(|V|V'fc|)) 

n— >oo ' 
A:=l 

liminf TV (|V|7„|V|). 
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We thus obtain the desired result. □ 

We are now in position to prove Theorem [D Let (7n)neN be a minimizing sequence for Ix- 
We know from Lemma [2] that (7n)neN is bounded in Ti and that {^p-y„)nm is bounded in 
H^{M.^). Replacing (7n)nGN by a suitable subsequence, we can assume that (7„) converges 
to some 7 S /C for the weak-* topology of TC and that (\/P7„)neN converges to y/p^ weakly 
in ff^(M^), strongly in L|'q^(M^) for all 2 < p < 6 and almost everywhere. 

If Tr (7) = A, then 7 G ICx and according to Lemma (HJ 

f^°^(7)<liminffLDA(^^)^j^ 



yielding that 7 is a minimizer of (|35p . 

The rest of the proof consists in rulling out the eventuality when Tr (7) < A. Let us 
therefore set a = Tr (7) and assume that < a < A. Following e.g. 0, we consider 
a quadratic partition of the unity + = where ^ is a smooth, radial function, 
nonincreasing in the radial direction, such that ^(0) = 1, < ^(x) < 1 if |x| > 0, ^(x) = 
if |x| > 1, IIV^IIl- < 2 and ||V(1 - < 2. We then set (,r{-) = For 

all n G N, 1-^ Tr (^Rjn^R) is a continuous nondecreasing function which vanishes at 
R = and converges to Tr (7„) = A when R goes to infinity. Let Rn > be such that 
^ {(.RnlnCRn) = Oi. The sequence {Rn)neN goes to infinity; otherwise, it would contain a 
subsequence {Rnk)ke'N converging to a finite value R*, and we would then get 

/ p-yix)fji*{x)dx = Urn / p^„(x)C|^ (x)dx = 2 lim Tr (^H„7nCfl„) = 2q = / p^{x)dx. 

As Cr* < 1 on \ {0}, 

we reach a contradiction. Consequently, (-R„)neN indeed goes to 
infinity. Let us now introduce 

7i,n = CR„7n^i?,„ and 72,n = XR„7nXi?„- 



< 



Note that 7i^„ and 72,n are trace-class self-adjoint operators on L (M ) such that < 7^-^ 
1, that p-y„ = p^^ ^ + /9^2 ^ and that Tr (7i,n) = ot while Tr (72,n) = A — a. Besides, using 
the IMS formula 



2 

n I ; 



-A = Xi?„(-A)x/?„ +Cft.(-A)Cfl„ - |Vxi?„P - \ViR 
it holds 

Tr (-A7„) = Tv (-A7i,„) + Tr (-A72,„) - Tr ((|Vxfij' + \'^iRf)ln) 

4A 

> Tr (-A7i,„)+Tr(-A72,„)--^, (49) 

from which we infer that both (7i,n)neN and (72,n)neN are bounded sequences of TL. As 
for all (/> G (:7c"^(R3), 

IV {liM)m) = Tr {ln{\iRA){iRM) 

= Tr (7n(|(ei?„ - l)0)(e«„<A|)) + Tr {^n{\<P)m„ - 1)01)) + Tr (7n(|</>)(</'|)) 
Tr(7(|<A)(<Al)), 

n— +00 

we obtain that (7i,n)neN converges to 7 for the weak-* topology of TL. Since Tr (7i,n) = 
a = Tr (7) for all n, we deduce from Lemma [6] that (/07i_„)nGN converges to strongly in 
LP(M3) for all 1 < p < 3, and that 

fLDA(^) < £:L°^(7i,n). (50) 
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As a by-product, we also obtain that (/072,„)neN converges strongly to zero in Lj'^^(]R^) for 
all 1 < p < 3 (since p^j,™ = P-iu - Pii,u with (p7„)neN and (P7i,„)neN both converging to 
in for all 1 < p < 3). Besides, using again (l49l) . it holds 

= Tr(-A7„)+/ p^^V + J{p^J+ [ g{p,J 

> Tr (-A7i,„) + Tr (-A72,„) + / P^,,„V + [ p^.^^V 

+ J (P7l,n ) + J (P72,„ ) + / 5(P7l,n + P72,n ] 



4A 



4A 

+ iiaip^i,^ + - 9{p-yi.„) - 9ip-y2,„)) - 



For R large enough, one has on the one hand 
and on the other hand 

{9{p'yi,„ + /372,n) - 9{P-Il.n) - 9{PJ2,J) 



I 



<2Zl 1^ p7,,„) ||V^/p^||i2 + 



2Z{X-a) 
R 



< 



+ 



I^W.n +P72,n) -5(/57l,n)| + / \9{p'y2.n)\ 

J Br 

I^W.n + P72,n) - 5(P72,n)| + / I^W.n)! 



< c 



(P72,n+P72,„) + ll^7l,nllL2 / /J^^.n 



+ ^(X/72-,„+/72t„) 



+ c 



W.n +P71„) + l|P72,nllL2 



/'7l,n 



+ / P'.l.+P'. 



Bi 



71, n 



for some constant C independent of R and n. Yet, we know that (^/'7i,„)neN and 
(■\//^Ti^)"6N are bounded in f/'-'^(M^), that (P7i_„)ngN converges to ^7 in ^^(M^) for all 
1 < p < 3 and that (P72,,i)ngN converges to in L[|^^(M^) for all 1 < p < 3. Consequently, 
there exists for all e > 0, some G N such that for all n > A^, 



Letting n go to infinity, e go to zero, and using (l37l) . we obtain that Ia = /q + /^^q and 
that (7i,n)nGN and (72,n)nGN are minimizing sequences for Iq and respectively. It also 
follows from ([501) that 7 is a minimizer for Iq,. In particular 7 satisfies the Euler equation 



7 = l{-oo,.p)(-f^pJ +'5 

for some Fermi level ep G M, where 



P7 



^^ + V + p^^\v\-^ +g'{p^), 
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and where < 5 < 1, Ran((^) C Ker{Hp^ — ep)- As V + p^-k |r| ^ + g'{p^) is A-compact, 
the essential spectrum of Hp^ is [0, +oo). Besides, Hp^ is bounded from below, 

Hp^ < -^A + y + p^*|rr\ 

and we know from [171 Lemma II. 1] that as — '}2ik=i + J^3 P7 = —Z + 2a < —Z + 2X < 0, 
the right hand side operator has infinitely many negative eigenvalues of finite multiplicities. 
Therefore, so has Hp^ . Eventually, ep < and 

n m 
i=l i=n+l 

where < nj < 1 and where 

-^A^i + V(l)i + {py -k |r|"^) (pi + g'{py) (pi = ei (pi 

ei < £2 < £3 < • • • < denoting the negative eigenvalues of Hp^ including multiplicities 
(by standard arguments the ground state eigenvalue of Hp^ is non-degenerate). It then 
follows from elementary elliptic regularity results that all the (pi, hence /j^, are in iJ^(M^) 
and therefore vanish at infinity. Using Lemma [121 the (pi decay exponentially fast to 
zero at infinity. 

Let us now analyze more in details the sequence (72,n)ngN- As it is a minimizing 
sequence for (/072,„)ngN cannot vanish, so that there exists r/ > 0, i? > and such 

for all n G N, Jy^^^ii P'y2,n ^ V for some y„ £ M^. Thus, the sequence {Ty,,'j2,nT-yJnen 
converges for the weak-* topology of Ti to some 7' G /C satisfying Tr (7') > > 0. Let 
/? = Tr (7'). Reasoning as above, one can easily check that 7' is a minimizer for and 
that Ix = Ia + + lyLa^fs- Besides, 

7' = l(-oo,4)(^py)+'^' 

where 

H^^,=-^A + pyk\r\-'+g'ipy), 

and where 0<6' <1, Ran{6') C Ker(if;^, - e^), and ep. < 0. 
Assume for a while that one can choose ep/ < 0. Then 

n' m' 
1=1 i=n'+l 



all the 0i's being in C°°(M ) and decaying exponentially fast at infinity. For n G N large 
enough, the operator 

7n = min (l, II7 + Tnej' T-ne\\'^) {l + TneV'^-ne) 

then is in /C and Tr (7„) < [a + 13). As both the (piS and the (p'^s decay exponentially 
fast to zero, a simple calculation shows that there exists some (5 > such that for n large 
enough 

n '"^ n 
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Hence, for n large enough 

W < < ^^°^(7n) < + 

Adding to both sides, we obtain that 

which obviously contradicts the previously established equality Ix = la + I"^ + IxLa-p- 

It remains to exclude the case when ep' has to be chosen equal to zero. In this case, 
is an eig envalue of H^, and there exists V S Ker{H^,) C H'^{R^) such that H-^IIl^ = 1 
and 'y'lp = flip with /x > 0. We then define for < ?] < /x and n € N, 

7„,,, = min (l, ||7 + r/|(/)m+i)((/)m+l| + T„e(7' - ??lV')(V'l)T-neir"^) 
(T + ??|0m+l)(0m+l| + Tneil' " ??| l/') (V'l )T-,^e) • 

As 'jn^T] is in IC and such that TV (7n,,j) < A, it holds 
It is then easy to show that 

n— >oo 

Besides, for r/ > small enough 

Reasoning as above, we obtain that for r/ > small enough 

Ix< h + 2r]em+i + o(??), 
which is in contradiction with the fact that 6^+1 is negative. The proof is complete. 

5.4 Proof of Theorem d 

For (j) £ i?i(R3), we set p^{x) = 2\(j){x)\^ and 

E{cP)= [ \V<P?+I p^V + J{p^) + Eg''^{p^). 
For all 4) G i7^(M3) such that ||(/>||i2 = 1,7^ = |(/))((/.| G /Ci and £:(7<^) = Therefore, 
h < inf e //^(M^), |<^|2 = i| . 

Conversely, for all 7 G /Ci, ^7 = y ^ satisfies c/)^ G //-"^(M^), ||(/>||i2 = 1 and 



^GGA(^) = + TV (-A7) - - / |V^p > £^''^{\4>,){<^,\) = E{ct>,) 



1 

2 

Consequently, 

h = inf G //^M^), |</>|2 = l| (51) 
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and (f20l) has a minimizer for A'p = 1, if and only if ([5T]l has a minimizer (j) (7^ then is a 
minimizer of (f20]l for Np = 1). We are therefore led to study the minimization problem 
([STIl . In the GGA setting we are interested in, i?((/>) can be rewritten as 

E{ct>)=f |V0|2+/ p^V + J{p^)+[ %^,|V0|2). 

Conditions ([29]l - ([33]) guarantee that E is Prechet differentiable on H^{W^) (see [T] for 
details) and that for all {<j),w) £ H\R^) x H\R'^), 

E'{ct>yw = 2Q (1 + 1^ (p^, |V0p)) v^.v^+J^^ (v + p^* |r|-i + 1^ (p^, |V0p)) . 
We now embed ([5T]) in the family of problems 

Ja = inf !^E{^), (t> G H\M.^), J^^ \4>\^ = a| (52) 
and introduce the problem at infinity 

Jr = inf|i?°°(,/.), ct>GH\R'), lj^\^ = xj (53) 

where 

E^{<P)=[ |V0|2 + J(p^)+ / M/'</.,|V0|2). 
JR3 Jrs 

Note that reasoning as above, one can see that Jx = Ix and = for all < A < 1 
(while these equalities do not a priori hold true for A > 1). 

The rest of this section consists in proving that ([52]) has a minimizer for all < A < 1. 
Let us start with a simple lemma. 



Lemma 7. Let < p < 1 and let ((/>n)n6N be a minimizing sequence for (resp. for J^) 
which converges to some cj) G //""^(M^) weakly in H^iR?). Assume that ||(/>||^2 = A*- Then (f) 
is a minimizer for (resp. for J^)- 



Proof. Let {4>n)neN be a minimizing sequence for which converges to (p weakly in 
For almost all x G R^, the function z 1— > \z\'^ + h{p^{x), |zp) is convex on M^. Besides the 
function 1 1-^ t+h{p^(x),t) is Lipschitz on M_|_, uniformly in x. It follows that the functional 



I (iwf + %0,|vv|')) 



is convex and continuous on H'^i^^). As {4>n)n<^n converges to (j) weakly in //-'^(M'^), we get 

/ (|V</<|2 + /i(p^,|V0|2)) < lim / (|V0„|2 + /i(p^,|V(/>„|2)). 
Besides, we deduce from (1311) that 



{h{p^^,\V(t)n?)-h{p^,\y4>n?)) 
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where the constant C only depends on h and on the bound of ((/>n)nGN- As {(t)n)n&i 
converges to (f) weakly in Lp'^E?) and as \\<P\\l'2 = \\(l)n\\L'^ for all n G N, the convergence of 
['Pn)n&^ to (j) holds strongly in L^(M^). Therefore, 

/ \v<t^\''+E^,^\p^) = [ {\v^\'+h{p^,m')) 

< hminf / (|V0„|2 + /i(p<^,|V0„|2)) 

+ lim / {hip^„,\Vcl)n\^)-hip^,\V(l)n\^)) 

= liminf / |V</>„|2 + i?GGA( j_ 

Finally, as {(j)n)nm is bounded in and converges strongly to (/> in L^(M'^), we infer that 
the convergence holds strongly in L^'(M^) for all 2 < p < 6, yielding 

Therefore, 

E{4>) < liminf =/^. 

n— >oo 

As ||</'||^2 is a minimizer for J^. Obviously, the same arguments can be applied to 

a minimizing sequence for JJf . □ 



In order to prove that the minimizing sequences for J\ (or at least some of them) are 
indeed precompact in L^(M^), we will use the concentration-compactness method due to 
P.-L. Lions [18]. Consider an Ekeland sequence {(j)n)n€N for (f52l) . that is a sequence 
{(j)n)n&^ such that 

VnGN, (pn&H^O^^) and / (t>l = X (54) 

Jr3 

lim E{4>n) = Jx (55) 
lim ^'(</.„)+e„</.„ = ini7-i(M3) (55) 

n— ^+00 



for some sequence (^n)ngN of real numbers. As on the one hand, \4>\ G H [M."^) and 
i?(|(/)|) = E{(j)) for all (/) G //^(M^), and as on the other hand, the function A Ja is 
decreasing on [0, 1] , we can assume that 

Vn G N, 0„ > a.e. on and 9^ > 0. (57) 

Lastly, up to extracting subsequences, there is no restriction in assuming the following 
convergences: 



weakly in H^(M.^), (5 

'P 

^loc 



strongly in Lf„„(K^) for all 2 < p < 6 (59) 



) a.e. in (60) 
in M, (61) 
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and it follows from ([571) that (p '>0 a.e. on and 6 > 0. Note that the Ekeland condition 
([56]) also reads 

-^div ((l + 1^ (/0<^n, IVcAnl^)) V0„,^ + (V + ★ |r|-l + 1^ |V</.„|2)^ 0„, + en<i)n 

(62) 



with r7„ — > in i7 



-1/Tn)3\ 



We can apply to the sequence (0n)nGN the following version of the concentration-compactness 
lemma. 

Lemma 8 (Concentration-compactness lemma |18j). Let A > and ((/>n)neN &e a hounded 
sequence in i7^(M^) such that 



Vn e N, 



0n = A. 



Then one can extract from {4'n)nef-i ^ subsequence {<pnk)k<^n such that one of the following 
three conditions holds true: 



1. (Compactness) There exists a sequence {yk)k&i in such that for all e > 0, there 
exists R > such that 

yk G N, 



Vk+BR 



2. (Vanishing) For all R > 0, 



lim sup 



k—>oo 



Jy+BR 



0. 



3. (Dichotomy) There exists < 6 < X, such that for all e > there exists 

• a sequence (yA;)fceN of points ofM.^, 

• a positive real number Ri and a sequence of positive real numbers {R2,k)keN 
converging to +oo, 

• two sequences {4'i,k)k£N o-nd {cl)2,k)neN bounded in H^{M.^) (uniformly in e) 
such that for all k: 

= 4>i,k on yk + Br^ 
4>nk = <p2,k on \ {Vk + Br^ J 



\k 



< e, 



<l>lk - (A - 5) 



< e 



lim dist(Supp 0i^k) Supp 02,k) = oo 



6— p 



I'/'nfe - {(t>i,k + (/'2,fc) ||lp{r3) < Cp e 2p for all 2<p <6 
MLPiy,+iBR^,^\BR^)) < Cpe^ for all 2<p<6 
liminf / (\V^nf - |V</.i,fe|2 - \V^2,k\^) > -Ce, 

where the constants C and Cp only depend on the bound of {(j)n)neN- 
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We then conclude using the following result. 



Lemma 9. Let {(j)n)nm satisfying [5^-^61]}. Then using the terminology introduced in the 
concentration-compactness Lemma\^ 

1. if some subsequence {(j)n^)k<m of {(j)n)nm satisfies the compactness condition, then 
(0nfc)fceN converges to (j) strongly in LP(M'^) for all 2 < p < 6 ; 

2. a subsequence of ((/'n)ngN cannot vanish ; 

3. a subsequence of (0n)neN cannot satisfy the dichotomy condition. 

Consequently, ((/>n)neN converges to strongly in LP(R^) for all 2 < p < 6. It follows that 
(j) is a minimizer to j'5^|) . 

As the explicit form of the functions (pi^k and (p2,k arising in Lemma [8] will be useful for 
proving the third assertion of Lemma [9j we briefly recall the proof of the former lemma. 



Sketch of the proof of Lemma El The argument is based on the analysis of Levy's concen- 
tration function 

Qn{R) = sup / (f)l. 

The sequence (Qn)nGN is a sequence of nondecreasing, nonnegative, uniformly bounded 
functions such that lim Qn{R) = A. 

There exists consequently a subsequence ((5nfc)fceN and a nondecreasing nonnegative func- 
tion Q such that {Qn^)k&i converges pointwise to Q. We obviously have 



lim Q{R) = 5 G [0,A]. 

R— ►oo 

The case 5 = d corresponds to vanishing, while 5 = \ corresponds to compactness. We 
now consider more in details the case when < 5 < A (dichotomy). Let ^, x be in C°°(M^) 
and such that < < 1, i{x) = 1 if \x\ < 1, i{x) = if |x| > 2, x{x) = if |x| < 1, 
X{x) = 1 if |x| > 2, IIVxIIl- < 2 and ||V^||loc < 2. For ii > 0, we denote by = ^ {^) 
and xr{') = X ("r)- Let e > and Ri > e^^ large enough for Q{Ri) > 5 — | to hold. 
Then, up to getting rid of the first terms of the sequence, we can assume that for all k, we 
have Qn^{Ri) > 6 — e and Qn^{2Ri) < (^ + |. Furthermore, there exists G such that 

Qn,{Rl)= j <tl, 

and we can choose a sequence {R'j^)k£N of positive real numbers greater than converging 
to infinity, such that Qn^,{2R'j^) < 5 + e for all A; G N. Consider now 



4>i,k = ^rA- - yk)(t>', 



and 



Denoting by i22,fc = we clearly have 



'}2,k 



blk 



< e, 



^Ik 



(A -5) 



< e, 



^Ik 



Ri<\--yk\<R2,k 



< QnkiR2,k) - QndRl) < 2e, 



26 



and 



/ 



^-i.fc + </'2,fc)|^ < / |i -yk)-XR'^{- -yk)?4>l^, 



< 



Ri<\-yk\<R2,k 



< 2e. 



Similarly, by Holder and Gagliardo-Nirenberg-Sobolev inequalities, we have for all k and 
2 <p < 6 

(6-p) 

^^< C„e 



+ (p2,k)\\LP < 



where the constant Cp only depends on p and on the bound on (0n)neN- Finally, we 
have llV^ijJIioo < 2/?j"^ < 2e and ||Vxr' < 2{R';^)-^ < 2e, so that 



|v</>i,fcp-eii(--yfc)|V0„J 



and 



N(p2,k\^ -x\{--yk)\v^nf 

where the constant C only depend on the bound on ((/)„)neN- Thus 

/ \V4>nf-\V4>i,k\'-\^<t>2,k\^ > I (i-eli(--yfc)-xk(--yfc))|v</.„j2_ce 



> -Ce. 



□ 



Proof of the first two assertions of Lemma\^ Assume that there exists a sequence {yk)ke'N 
in M.^, such that for all e > 0, there exists R > such that 

yken, [ > A - e. 

Two situations may be encountered: either {yk)keN has a converging subsequence, or 
lim l^fcl = cxD. In the latter case, we would have = 0, and therefore 

fc^oo 



Hence 



lim 

fc— >oo 



< lim E^icP^^: 

fe— >oo 



V = 0. 



lim E{ 

k^oo 



which is in contradiction with the first assertion of Lemma ([T]). Therefore, {yk)keN has a 
converging subsequence. It is then easy to see, using the strong convergence of {(j)n)neN to 
that 



> 



> A 



V+Br 



where y is the limit of some converging subsequence of {yk)keN- This implies that \\4>\\\2 
A, hence that {4'n)neN converges to (j) strongly in L^(M^). As ((/>n)ngN is bounded in 
this convergence holds strongly in LP(M^) for all 2 < p < 6. 

Assume now that {(j)nk)k£n is vanishing. Then we would have (/> = 0, an eventuality that 
has already been excluded. □ 
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Proof of the third assertion of Lemma\^ Replacing ((/>„)neN with a subsequence and using 
a diagonal extraction argument, we can assume that in addition to (f54ll - ([6T]) . there exists 

• a sequence {yn)n&i of points in M^, 

• two increasing sequences of positive real numbers (i?i,n)nGN and (i?2,n)neN such that 

lim Rin = oo and lim R2n — Rin = oo 

• two sequences (</>i,n)n6N and (02,n)neN bounded in i7^(M^) 
such that 

i>n = (l>i,n on y„ + „ 

^n = 02,n on \ {Vn + Br^ J 



lim 

n— >oo 



lim 

n— >oo 



lim ||0„ - (01 „ + (f)2,n)\\Lp(m = for all 2 < p < 6 

I— >oo ^ ' 

\\Lr'(y„+iBn,,^BR,,J) =^ for all 2 < p < 6 



lim 

n^oo 



lim dist(Supp (pi n, Supp i?i'2 n) = co 
n^oo ' ' 

liminf / (\V<Pnf - \y<Pl,nf - |V(/<2,n|') > 0. 

Besides, it follows from the construction of the functions (pi^n and (/)2,n that 

Vn e N, > and (/)2,n > a.e. on M^. 

A straightforward calculation leads to 



(63) 



+ 



(|V0„|2-|V0i,„p-|V</.2,„|2) + / A 

3 V / ./HPS 



y 



+ [ (^(/'0„,|V</.„|2)-/l(p<^^_„,|V0i,„|2)-/i(p^ |V02,n|')), (64) 



where we have denoted by pn = p 



As 



the sequence (pn)n6N goes to zero in LP(R'^) for all 1 < p < 3, yielding 



R3 



ni Pr, 



0. 



Besides, 



D{P4>i,„,P<j>2,u) < 4 dist(Supp (/)i,n,Supp (^2,n) Ml</'l,n|li2 ||</'2,n|li2 > 
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and 



< 



|V0„|2) _ |V(/.i,„P) - |V02,nP)) 



< c 



(recall that 1 < p± = 1 + /9± < |). Lastly, as Iim„^oo dist(Supp ^i^n, Supp (/)2,n) = oo, 



mm 



0. 



It therefore follows from (l64l) and from the continuity of the functions A i-^ Ja and A i-^ J' 
that at least one of the inequalities below holds true 



oo 
A 



J\> Js + JT~5 (case 1) or J\> Jf" + Jx^s (case 2). 
As the opposite inequalities are always satisfied, we obtain 

Jx = Js + Jx^s (case 1) or Jx = + Jx-5 (case 2) 
and that (still up to extraction) 



(65) 



(66) 



lim E{(l)in) = Js 

lim i?~(02,n) = 



hm 



(case 1) or 



lim E((j)2,n) = Jx-5 



(case 2). (67) 



Let us now prove that the sequence ('0n)neN, where i^n = (t^n — {(t^i,n + </'2,n), goes to zero 
in H^{R^). For convenience, we rewrite ipn as ipn = encf^n where Cn = 1 — C-Ri,„(" ~ Vn) — 
X_R2n/2(" ~ Vn) and Ekeland's condition ([62]) as 



div {anV(t^n) + V(t>n + (/>^„ * \A~^)(t>n + KT^^'^" + + ^nK = Vn (68) 



where 



1 / 5/i , 1^ , ,n. 



,>i- 



K+ = 2^V;^^(p0„,|V0„|^)Xp,„ 

The sequence (y(/)„ + {p^^ -k \r\~^)(f>n + V~(j)]^'^^' + V^<j)]^'^'^* + 6'„(?!)„)„gpj is bounded in 
L^(M^), (r7„),igN goes to zero in H~^{E?), and the sequence ('i/'n)ngN is bounded in H'^i^^) 
and goes to zero in L^(M^). We therefore infer from (f68ll that 



n— »oo 



Besides VV'n = e„V(/)„ + (/)„Ve„ with < < 1 and ||Ve„||Loo — > 0. Thus 



0. 
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As 



a if dh , ,„ .nA h „o , , 

0<77<an = 7T l + 7rW.JW„|2) <-<oo a.e. onM^ (gg) 



and < < e„ < 1, we finally obtain 



/ el\V<t)n\'' 0, 

Jr3 n^oo 



from which we conclude that (VV'n)neN goes to zero in f/'^(M^). Plugging this information 
in ([68]) and using the fact that the supports of 4)i^n and (/>2,n are disjoint and go far apart 
when n goes to infinity, we obtain 

-div (a„V(/.i,„) + + {p^, „ * \r\-^)Kn + ^."0} n'^~ + K^'^lt'^-' + ^n-/)!,™ ^ 

' ' ' n— >oo 

-div (a„V02,n) + V^</'2,n + * \A'^)<t>2,n + Vn (t\n^~ + Vrt ^^n^^ + Gn^2,n 

We can now assume that the sequences {(j)i^n)n£N and (02,n)neN, which are bounded in 
H^{M.^), respectively converge to ui and U2 weakly in /f^(M^), strongly in LJ'q^(]R'^) for all 

2 < p < 6 and a.e. in M^. In virtue of (f63]l . we also have ui > and U2 > a.e. on 
M^. To pass to the limit in the above equations, we use a H-convergence result proved in 
Appendix (Lemma [TO]) . The sequence (a„)„gN satisfying (f69]) . there exists Coo S -L°°(R^) 
such that I < Ooo < |^ and (up to extraction) a„/3 -^h ciooH (where I3 is the rank- 

3 identity matrix). Besides, the sequence {V^)n&i is bounded in L°°(R^), so that there 
exists G L°°(R^), such that (up to extraction) {V^)n<m converges to for the weak-* 



topology of L°°(M ). Hence for j = 1, 2 (and up to extraction) 
' V(l)in — > Vuj strongly in H'^m.^) 

n^oo 

Vn^cP]!''- - V^u]^'^- weakly in " 



j,n n^oo 3 weaKiy in ^1^^,^ 

3,n 



*|r| )(t)j,n + On(i)j,n ^ * |r| )uj + 9uj stronly in Lio^(M-'). 

n — >oo 

We end up with 

-div (oooVui) + Vui + {pu, ★ |rr^)Mi + V'ul'^^'^- + y+n}+^^+ + Qui = (70) 
-div (oooVua) + Vu2 + ipu2 * \r\'^)u2 + V'ul^^'^- + V+ul'^^^+ + = 0. (71) 

By classical elliptic regularity arguments [Qj (see also the proof of Lemma [T2l below) . both 
ui and U2 are in C'''"(M^) for some < a < 1 and vanish at infinity. Besides, exactly one 
of the two functions ui and U2 is different from zero. Indeed, if both ui and U2 were equal 
to zero, then we would have (j) = 0, hence 



TOO 



Ja = lim E{<Pn) = lim (</<„) = , 

which is in contradiction with the first assertion of Lemma [U (recall that Jx = I\ and 
= for all < A < 1). On the other hand, as dist(Supp 0i^n;Supp </>2,n) — ^ oo, at 
least one of the functions ui and U2 is equal to zero. 

We only consider here the case when U2 = 0, corresponding to case 1 in (f65l) - (f67l) . 
since the other case can be dealt with the same arguments. A key point of the proof 
consists in noticing that apply Lemma [TT] (proved in Appendix) to ([TOl) (note that W = 
V'uf'^ + V'^u^'^ is nonpositive and goes to zero at infinity) yields 

^ > 0. (72) 
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Consider now the sequence {(pi,n)neN defined by (j)i^n = '^^</'i,n||</'i,n||^2^- It is easy to check 
that 

Vn G N, ^1 „ G H\M.^), / „ = 5 and n > a.e. on 
lim E{(j)i^n) = Js 

-div (ai,„V^i,„) + + (pT * \r\-^)^l,n + ^fn^lt'^" + ^I'^n^lt'^^ + ^n^l.n ^ 

{4'i,n)neN converges to 7^ weakly in H^, strongly in L^^^ for 2 < p < 6 and a.e. on M 

(with in fact vi = (/>). Likewise, the sequence ||</'2,n||^2 4'2,n)n€N being a minimizing 

sequence for J^^, it cannot vanish. Therefore, there exists 7 > 0, i? > and a sequence 

{xn)nm of points of such that \4'2,n\'^ > 7- Then, denoting by (^2,71 = (A — 

(5)2 ||02,n||^2 4'2,ni- " 

VnGN, 4>2n^H\M.^), / 0^ ^ = A - 5 and ^2 n > a.e. on 

hin i?-(02,n) = Ja°°_5 
< n— »+oo 

-div (a2,„V^2,n) + * |rr^)02,n + ^2":n?2r''" + ^2^^'^^ + ^"^2,n^^0 

('?^2,n)neN Converges to U2 7^ weakly in H^, strongly in Lf^^ for 2 < p < 6 and a.e. on R 
It is important to note that the sequence {aj^n)nm and (^^)neN are such that 

^ < aj,n < ^ and < 2^+C, 

where the constants a, 6 and C are those arising in (f3T]l and ([33il . 

We can now apply the concentration-compactness lemma to {<f>i^n)neN and to ((/'2,n)neN- 
As ((/>j,n)neN does not vanish, either it is compact or it splits into subsequences that are 
either compact or split, and so on. The next step consists in showing that this process 
necessarily terminates after a finite number of iterations. By contradiction, assume that 
it is not the case. We could then construct by repeated application^ of the concentration- 
compactness lemma (see [1] for details) an infinity of sequences {ipk,n)n€N, such that for 
all A; e N 

Vn G N, V'fc n e i?^(R^), / i^ln = ^k and V'fc n > a.e. on R^ 

{tpk,n)nm converges to Wk weakly in H^^ strongly in L^^^ for 2 < p < 6 and a.e. on 1 
with 

Y,6k<X, (73) 

fceN 

and with for all A: G N, 

I < Sfc,„ < ^ and ||i^fc=;=„||L- < 2^+C. 

Using Lemma fTOl to pass to the limit with respect to n in the equation satisfied by V'fc,n, 
we obtain 

- div {akVwk) + {p^, * \r\-^)wk + V^w^r""- + K^l"'^"^ + Owk = 0, (74) 
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with 

a ~ 
-<Sfc<— and WV^Wl^ <2f^+C. 

Besides, we infer from ((73]) that ||'W^fc|li2 < A, hence that 

km 

Hm ||ii'fc||L2 = 0. 

It then easily follows from ([74| that 

lim ||div {akVwk)\\L-i = 0. 

fe— >oo 

We can now make use of the elliptic regularity result [9] (see also the proof of Lemma [T2ll 
stating that there exists a constant C, depending only on the positive constants a and h, 
such that for all u G H^{M?) such that div (ofcVn) G L'^{M?), u G L°°(M3) and 

\\u\\l^ < C{\\u\\l2 + ||div {akVu)\\L2) 

and obtain 

lim ||u)fc||ioo = 0. 
k^oo 

Lastly, we deduce from (f74l) that 

n\\ l|2 / ^ / II ||2/3- 2/3+\ ,, ,,2 

0\\wk\\L2 < C [\\wk\\£oo + \\wk\\£oi j ||^^fc|li2. 
As ||i(;fc||i2 > for all fc G N, we obtain that 

9 <C (\\wk\fj^c^ + lliffclli^ ) — 0, 

\ / k—*oc 

which obviously contradicts (f72l) . We therefore conclude from this analysis that, if di- 
chotomy occurs, {(l)n)neN splits in a finite number of compact bits. We are now going to 
prove that this cannot be. 

If this was the case, there would exist 6i > and ^2 > such that < 6i + 82 < X and 
two sequences (iii,n)nGN and (u2,n)neN such that 

VnGN, ni,„ G i^^(M^), / jiii,™!^ = ui > a.e. on 
< lim E{ui^n) = hi 

n— >oo 

-div (ai,„Vni „) + Vui ^ + (Pm „ * |rr"^)^il,n + t^fn^lt^'^" + ^^n^lt^'^^ + ^nUl,n ^ 

and 

VnGN, U2„Gi/^(M^), / |n2„|^ = 52, ii2 > a.e. on 

< lim E'^{u2,n)=l52 
n—>oo 

-div (Q2,nVn2,„) + {pu2„ -k \r\-'^)u2,n + V2„ul't^'^- + V^„U2t,^^+ + 6l„U2,n ^ 

and converging weakly in to ui and U2 respectively, with ||mi||l2 = 61 and ||w2||l2 = 

82 (as the weak limit of {(j)n)neN in L^(M^) is nonzero, one bit stays at finite distance from 
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the nuclei). It then follows from Lemma [7] that Ui and U2 are minimizers for J^^ and 
respectively: 

E{ui) = Jsi, ||tii||^2 = 5i, E{u2) = J^, \\u2\\l2 = 62- 
Letting n go to infinity in the equations satisfied by ui^n and U2^n we also have 

- div (oiVni) + Vui + {pu, ★ \r\^^)ui + v^ul-^^^' + vtul'^^^+ + Oui = (75) 

and 

— div (a2Vn2) + (puj * |r|^"'^)u2 + 1^2^ Ug'*'^^" + ^^1*2^^^+ + 9u2 = 0, (76) 

with ^ < aj < ^ and ||t;j'^||ioo < This shows in particular that ui and U2 are 

in L°°(R^). Applying Lemma [T2l we then obtain that there exists 7 > 0, /i € //^(M^), 
/2 e i^HlR^), 91 e (^^(M^))^ and 52 e (^^(M^))^ such that 

Ml = e~'^'''/i, 'U2 = e"'^'''/2, Vtii = e"'^l'ls'i, Vu2 = e"'^l'l5'2- (77) 

In addition, as ui > and U2 > 0, we also have /i > and /2 > 0. Let e be a given unit 
vector of M^. For i > 0, we set 

wt{r) = at {ui{r) + U2{r - te)) where a* = ((5i + (^2) 2 + U2(- - te)||^2 . 

Obviously, G H^{R^) and ||tfi||L2 = + ^2, so that 

E{wt) > Js,+5,. (78) 

Besides, 

\\ui+U2{--te)\\l, = [ ul+ [ 4 + 2 [ /i(r)/2(r-te)e-^(lrWr-*-l)dr 

JR^ JR3 JIR3 

= Si + 62 + 2 [ /i(r)/2(r-te)e-^(H+|r-te|)^r 

JM3 



yielding 

Likewise, we have 



di + 62 + 0{e-^'), 

at = l + 0(e-^*) 



/ \Vwt\^= [ |Vni|2+ / |Vt.2|' + 0(e-^*) (79) 

JR3 JR3 JR3 

/ V\wtf=[ V\ni\^+[ V\u2{--te)\' + 0{e-^') (80) 

Jr3 Jm3 Jir3 

D{pwt,Pwt) = D{pui,pu^) +D{pu2,Pu2) + '^D{pu^,p^^(^._te)) +0(e"^*). (81) 
The exchange-correlation term can then be dealt with as follows. Denoting by 

n = Pwt - Pui - Pu2{-te) = 2(at - l)(hi|^ + \u2{- - te)f) + AatUiU2(: - te) 

and 

St = \Vwt\^-\Vuif-\Vu2{--te)f = {af-l){\Vuif + \Vu2{--te)\^)+2a^tVuvVu2{-te), 
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and using ([3T]) . ([331) . ([771) and the fact that ui and U2 are bounded in we obtain 

HPrnA^wtl"^) - |Vni|2) - /i(p„2(._te), iVuzl- - te)|2) 

< / |/i(/0«i +Pn2(-te) |Vui|^ + |Vu2(- - + Sj) - |Vni|2)| 

2 

+ / |^(P«2(-te) + P«i + n, |Vn2(- - te)|2 + |Vui|^ + St) - h{pu2{-te),\^U2{- - te)\^) 

Jte+Bt 

+ / |/i(/>,2(.-te),|VtX2(--te)|2)| + / |/l(p,,,|V7/i|2)| 
JBt Jte+Bt 



+ I / |M/5^,,|Vu;t|^)|+/i(p,„|Vtxi|^)| + |%,,(._te),|Vn2(--te)|^)| =0(e-^*) 



BtyJ{te+Bt ) 
1. 1. 



Combining (I79|1 - (|8T]) together with the above inequahty, we obtain 

E{wt)<J5, + Jr, + I V\u2{--te)\^ + D{pu,,Pu,i.-te)) + 0{e-^'). 
Next, using (f77l) . we get 

/ yPu2{-te)+ D{pu^,Pu,^(^._te)) = -Zt~^ Pu2+t~^ Pui / p«2 + 

= -2J2(^-2Ji)t-^ + o(t~^). 

Finally, 

for t large enough, which contradicts (fTSl) . □ 



£?nd o/ ^/le proof of Lemma O As a consequence of the concentration-compactness lemma 
and of the first three assertions of Lemma [H the sequence {(j)n)n€N converges to (p weakly 
in i?^(R^) and strongly in LP(R^) for all 2 < p < 6. In particular, 

- hm / (Pl = \. 



n— >oo 



It follows from Lemma [7] that cj) is a minimizer to (f52l) . □ 

Appendix 

In this appendix, we prove three technical lemmas, which we make use of in the proof of 
Theorem [2l These lemmas are concerned with second-order elliptic operators of the form 
— div (^V-). For the sake of generality, we deal with the case when A is a matrix- valued 
function, although A is a real-valued function in the two-electron GGA model. 

For O an open subset of and < A < A < oo, we denote by M{X, A, J7) the closed 
convex subset of L°°(r2,M^^^) consisting of the matrix fields A G L°°(r2,M^^^) such that 
for all S and almost all x G 0, 

A|eP<^(xK-C and \Aix)C\<MC\. 
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We also introduce the set M^(A, A, Q) of the matrix fields A € M{\, A, fi) such that A{x) 
is symmetric for almost all x G 17. Obviously, M'^(A, A, J7) also is a closed convex subset 
of L°°(17,M3><3). 

The first lemma is a H-convergence result, in the same line as those proved in the 
original article by Murat and Tartar [T^, which allows to pass to the limit in the Ekeland 
condition ([62]) . Recall that a sequence (^n)neN of elements of M(A,A, 17) is said to H- 
converge to some A S M{X', A', Q), which is denoted by An A, if for every u CC 17 the 
following property holds : V/ G H~^{u}), the sequence (u„)„gn of the elements of Hq{lu) 
such that 

-div(^„Vn„) = fl in H-\uj) 

satisfies 

Un u weakly in Hq{uj) 
AnVun AVu weakly in L'^{uj) 

where u is the solution in Hq{lo) to — div(^Vu) = f\^. It is known that from 
any bounded sequence {An)n<^N in M(A,A, 17) (resp. in M(A, A, 17)) one can extract 
a subsequence which H-converges to some A G M(A, A~^A^, 17) (resp. to some A G 
M*(A,A-1A2,17)). 

Lemma 10. Let 17 be an open subset o/M^, < A < A < oo, < A' < A' < oo, and 
{An)n&N 0, sequence of elements of M{\, K,VL) which H-converges to some A G M(A', A',17). 
Let {un)n&%, {fn)neN dnd {gn)n€N be sequences of elements of H^{Q), H~^{Q) and L^(17) 
respectively, and u G H'^{^), f G //~"'^(17) and g G L^(17) such that 

— div(y4„Vn„) = fn + 9n in H^^{0,) for all n £N 
Un u weakly in H^{Q) 
fn^f strongly in H (17) 
^ 9n ^ 9 weakly in L^(17). 

Then — div {AVu) = f + 9 and AnVun AVu weakly in L^(17). 

The second lemma is an extension of [171 Lemma ILl] and of a classical result on the 
ground state of Schrodinger operators [26]. Recall that 

L'^{R^)+L'^{R^) = |w|Ve > 0, 3(W2,Woo) G l2(m3) X L~(r3) s.t. 




Lemma 11. ie^ < A < A < oo, A G M'{X,A,R'^), W G L'^{R^) + Lf=(M3) such that 
W+ = max(0, VF) G L^(M^) + L'^(]R'^) and fx a positive Radon measure on R^ such that 
/x(M3) < Z = Efeli^fc- Then, 

H = -div {AV-) + V + ni^\r\-^ + W 

defines a self-adjoint operator on L^(M'^) with domain 

D{H) = {n G F^(M^) I div (AVu) G l2(m3)} . 

Besides, D{H) is dense in H^{R^) and included in L°°(M3) n C°'"(M3) for some a > 0, 
and any function of D{H) vanishes at infinity. In addition. 
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1. H is bounded from below, o"ess(-ff) C [0, cxd) and H has an infinite number of negative 
eigenvalues; 

2. the lowest eigenvalue fxi of H is simple and there exists an eigenvector ui G D{H) 
of H associated with fii such that ui > on M^; 

3. if w £ D{H) is an eigenvector of H such that w > on M^, then there exists a > 
such that w = aui. 



The third lemma is used to prove that the ground state density of the GGA Kohn-Sham 
model exhibits exponential decay at infinity (at least for the two electron model considered 
in this article). 

6 

Lemma 12. Lei < A < A < oo, ^ e M{\,k,W^), V a function of L^^^iM?) which 
vanishes at infinity, 6 > and u G Lf^(M'^) such that 

-dW{AVu)+Vu + eu = inV'{R^). 

Then there exists 7 > depending on {X,A,9) such that e^'^^^u G H^{M.^). 



Proof of Lemma\T^ Let us denote by C„ = AnVun- One can extract from the sequence 
(Cn)neN, which is bounded in L^, a subsequence {£,nk)k£n which converges weakly in Lp'{Q) 
to some C solution to — div (C) = / + 5 in Lf~^(f2). The proof will be completed if we can 
show that we necessarily have ^ = AVu. Consider w CC O, g G H~'^{u)) and Vn G HQ{ijj) 
satisfying 

-d\Y{A*nVvn) =q in H'^{uj). 
As the sequence {A^)n£n //-converges to A* [19], it holds 

j V in Hq{uj) 

\ AlVvn A*Vv in L'^{uj) 

where v is the solution to —div {A*\/v) = q in Hq{lo). Let (p G C^{lo). As 

(f)Vn (t>v in Hq{uj) 
4>Vn — > i'V in LF'{uj) 
V(f)Vn V(i)v in {L^{uj)f 
^VcpUn^ Vcpuin {L'^{uj)f, 

we have on the one hand 



-(div^nk,</'^'nJ/^-i(a;),//iH - / ^"fe " V(/) 

J UJ 

J UJ J UJ 

-{fAv)H-Huj),H^{co)- / 9(t>V- I ^•V(/>V 

J UJ J UJ 

-(div^,(/.r;)^-i(^) - / ^■Vcl)v= / ^-Vvcf), 
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and on the other hand 



= - UnM*VVn^) • V0+ / U. 

J UJ J LU 



Therefore, 



/ u{A*Vv) • V(/. + / uq^= I Vn • {A*Vv)(^ = f (AVu) ■ V 
/ ^ . Vv0 = / {AVu) ■ Vv(f). 

J UJ J UJ 



As the above equahty holds true for all all v G Hq^lo) and all (p € C^(tj), we finally 
obtain ^ = AVu. □ 



Proof of Lemma \Tli The quadratic form go on L^(M'^) with domain D^qo) = 
defined by 

V(u,?;) G Z)(go) X ^(go), qo{u,v) = AVu-Vv, 



is symmetric and positive. It is also closed since the norm y\\ ■ ||^2 +'?o(") is equivalent 

to the usual norm. This implies that go is the quadratic form of a unique self-adjoint 
operator Hq on L^(M^), whose domain D{Hq) is dense in H^(M.^). It is easy to check that 
D{Ho) = {ue H^{R^) I div (AVu) E L'^{R^)} and that 

Vu e D{Ho), Hqu = -div {AVu). 

Using classical elliptic regularity results [9], we obtain that there exists two constants 
< a < 1 and C G M+ (depending on A and A) such that for all regular bounded domain 
17 CC M^ and all v G H^{Q) such that div (AVv) G L'^(n), 

Ibllco.-(n) := sup|u| + sup — ■ -— — < C(||?;||i2(f^) + ||div (^Vv)||i2(f^)) . 

n (r,r')ef^xQ l^^ " I 

It follows that on the one hand, D{Ho) ^ L°°{M.^) n C°'"(M3), with 

yueD{Ho), ||n||ioo(R3) + sup kW-yi <c(||^||^, + ||ffp^||^,), (82) 

and that on the other hand, any u £ D{Hq) vanishes at infinity. 

Let us now prove that the multiplication byW = y-|-^*|r|^^-|-M/^ defines a compact 
perturbation of f^o- For this purpose, we consider a sequence {un)n&i, of elements of D[Hq) 
bounded for the norm || • \\hq = (|| • ||^2 + ll-f^o " llia)^ . Up to extracting a subsequence, we 
can assume without loss of generality that there exists u G D{Hq) such that: 

Win ^^(M^) and LP{W^) for 2 < p < 6 
n„ ^ n in Lf^^(M^) with 2 < p < 6 
n„ — > n a.e. 



37 



Besides, it is then easy to check that the potential W = V + fi * \r\ ^ + W belongs to 
L2 + Lf (M^). Let e > and (W2, Woo) G L'^{M.^) x L°°(R3) such that ||Woo||l- < e and 
W = W2 + Woo. On the one hand, 

\\W^{un - u)\\l2 < 2e sup ||u„||ho, 

raeN 

and on the other hand 

lim \\W2{Un-u)\\L2=0. 

n— >oo 

The latter result is obtained from Lebesgue's dominated convergence theorem, using the 
fact that it follows from (f82l) that {un)nm is bounded in Consequently, 

lim \\Wun-Wu\\L2 = 0, 

n— >oo 

which proves that W is a ffo-compact operator. We can therefore deduce from Weyl's 
theorem that H = Hq + W defines a self-adjoint operator on L^(R^) with domain D{H) = 
D{Ho), and that acss{H) = aess{Ho). As qo is positive, (t{Hq) C M+ and therefore 

Let us now prove that H has an infinite number of negative eigenvalues which forms an 
increasing sequence converging to zero. First, H is bounded below since for all v G D{H) 
such that \\v\\i2 = 1, 



{v\H\v) = [ ^Vt;-Vt;+ / Wt;^ 

> X\\Vv\\l2 - \\W2\\l4'^v\\12 - e 



> -^A-3||w.ir 



In order to prove that H has at least negative eigenvalues, including multiplicities, we 
can proceed as in the proof of [17| Lemma ILl]. Let us indeed consider radial functions 
(pN in P(M^) such that for all 1 < i < A^, supp((/)j) G \ Bi and f^s = 1. 

Denoting by (j)i,a{-) = we have 



Jr3 



and 



1 1 

2 / p \ 2 
'2^11 T;r/2 \ II i. II 2 I ( / TT/3 1 II i. ||2 



1_ 

2 / ^ \ 2 



= o{a) 

where we have split W+ = max(0, VF) as W+ = W2 + W3 with W2 G L'^{M.^) and W3 G 
l'^). Besides, we deduce from Gauss theorem that 

Jr3 ' Jr^ max(|r|, |r'|) 7]r3 |r| 
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and that, for a small enough, 



Ml 2 



Thus, 



{(t)^,„\H\(t)i,„) < a{^{W^)-Z) [ M£^dr+ o (a), 

7r3 |r| 

yielding (^j^o-l-f^l^i.o-) < for cr > small enough. As (pi^a- and (pj^a have disjoint supports 
when i 7^ j , we also have 

max {4>\H\4>) < 

0espan((/)i,„,- ,9i]v,<T), Il</'lli2=i 

for C7 > small enough. It follows from Courant-Fischer formula [26] and from the fact 
that (Tess(-ff) C M+ that -ff has at least N negative eigenvalues, including multiplicites. 

The lowest eigenvalue of which we denote by fii, is characterized by 

/ii=inf|/ AVu-Vu+ W|np, ueH^R^), ||n||i2 = ll, (83) 



and the minimizers of (|83]) are exactly the set of the normalized eigenvectors of H associated 
with /ii. Let -ui be a minimizer (l83|). As for all u £ H^(M.^), \u\ G H^{R^) and V|n| = 
sgn(u)Vn a.e. on M^, \ui\ also is a minimizer to ((83]) . Up to replacing ui with there 
is therefore no restriction in assuming that ui > on M'^. We thus have 

ui G i?^(M^) nC°(M^), ni>0 and - div (^V-ui) + c/ui = 

with g = W — jJi G L^^^(M^) for some p > ^ (take p = 2). A Harnack-type inequality due 
to Stampacchia [28] then implies that if ui has a zero in M^, then ui is identically zero. 
As 1 1 ui 11^2 = 1, we therefore have ui > on M^. 

Consider now w G D{H) \ {0} such that Hw = fiw and w > on M.^. It holds 
/i / wiw = {w\H\wi) = /Ui / wiw. 

As w is not identically equal to zero and as ii^i > on M'^, J^^^wiw > 0, from which 
we deduce that fi = iii. It remains to prove that /ii is a non-degenerate eigenvalue. By 
contradiction, let us assume that there exists v S D{H) such that Hv = iiiv, \\v\\i2 = 1 
and {v,ui)i^2 = 0. Reasoning as above, also is an eigenvector of H associated with 
and \v\ > on M^. Since D{H) C (7'^(IR^), v is continuous on M^, so that either v = \v\ 
on or = —\v\ on M^. In any case, | Jjgs Uiv\ = Jj^s ui\v\ > 0, which is in contradiction 
with the fact that {ui,v)i2 = 0. The proof is complete. □ 

Proof of LemmaUM. Consider R > large enough to ensure 

on 

-<V{T) + e<— a.e. on5|j:=M3\5^. 

It is straightforward to see that u is the unique solution in H^{B'j^) to the elliptic boundary 
problem 

J -div(AV^;) + Vv + ev = m B% 

1 f = u on dBji. 
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Let 7 > 0, u = nexp "^^I'l and w = u — u. The function w is in ff^(]R^) and is the 
unique solution in H^{Bfj) to 



j -div(^V'u;) + + 6it(; = div(^V-u) -Vu-Ou in 5^ 
1 w = on OBr. 

Let us now introduce the weighted Sobolev space Wq{B'^^) defined by 

endowed with the inner product 

= / ^ e'^l'''(i;(r)u;(r) + Vu(r) • Vu;(r)) dr. 

Multiplying ([84l) by i?!)e^^l ' with (/> G T>{B'j^) and integrating by parts, we obtain 



(84) 



and then 



27|r| 



Jb^ Jb^ |r| Js^ 



(85) 



Due to the definitions of Wq{B'^) and u, ([85]l actually holds for (ti;,(/)) G Wo(-B^) x 
14^0^(5^), and it is straightforward to see that (f85ll is a variational formulation equivalent 
to (IMD . 

It is also easy to check that the right-hand-side in ((851) is a continuous form on WQ^(i?^), 
so that we only have to prove the coercivity of the bilinear form in the left-hand-side of 
([85]) to be able to apply Lax-Milgram lemma. We have for v G WqIBj^) 

[ ^e^l'"IVt;-e^l'"IV?; + 27 /" Ae^^^'^Vv ■ ^e^^^'^v + [ (V + e)e^ 
JB-h Jb"^ |r| 



> A 



2A7 



Bi 



+ 



> A 



>(A-A7) 



,7|r|v, 



2A7 



+ o 
L2{iJ^) '2 



A7) 



LHB-) ^ 2 



2 



2 



Thus the bilinear form is clearly coercive if 7 < min(-^, ^), and there is a unique w solution 
of in Wq{B'^) for such a 7. Now since u = w -|- u, it is clear that e'^' ln G H^{B'^), 
and then eTl l« G iJifM^). □ 
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